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We propose a microscopic theory of superconductiv¬ 
ity for systems with strong electron correlations such 
as cuprates in the framework of the extended Hub¬ 
bard model where the intersite Coulomb repulsion and 
electron-phonon interaction are taken into account. The 
Dyson equation for the normal and pair Green functions 
for the Hubbard operators (HOs) is derived. Due to 
the unconventional commutation relations for the HOs, 
a specific kinematical interaction of electrons with spin 
and charge fluctuations with a large coupling constant of 
the order of the kinetic energy of electrons W emerges 
that results in the d-wave pairing with high-Tc. Super¬ 
conductivity can be suppressed only for a large intersite 
Coulomb repulsion V > W. Isotope effect on Tc caused 
by electron-phonon interaction is weak at optimal doping 
and increases at low doping. The kinematical interaction 
is absent in the spin-fermion models and is lost in the 
slave-boson (-fermion) models treated in the mean-field 
approximation. 


INTRODUCTION 


To explain unconventional properties of cuprates one 
should take into account that cuprates are the Mott- 
Hubbard (more accurately, charge-transfer) doped insu¬ 
lators which cannot be described within the conventional 
band theory (for a review see, e.g., fl)- Under doping, 
a two-subband strongly-correlated metal emerges where 
the Fermi-liquid model fails to describe electronic excita¬ 
tions. The projected-type (Hubbard) operators referring 
to the two subbands must be introduced. A new energy 
scale of the order of the kinetic energy of electrons W 
arises in the intraband hopping induced by the kinemati¬ 
cal interaction for the HOs which is much larger than the 
antiferromagnetic exchange interaction J induced by the 
interband hopping proposed by Anderson Q. As shown 
in recent experiments ii, short-range antiferromag¬ 
netic (AF) dynamical spin fluctuations survive in super¬ 
conducting state even in the overdoped compounds. This 
justifies the spin-fluctuation mechanism of superconduc¬ 
tivity proposed earlier within spin-fermion models (see, 
e.g., Refs. d-Q). We consider the spin-fluctuation mech¬ 
anism of pairing induced by the kinematical interaction 
where the coupling constant is given by hopping param¬ 
eters. 


KINEMATICAL INTERACTION IN THE 
HUBBARD MODEL 

To describe electronic systems with strong correlations 
the Hubbard model is commonly used We consider 
the extended Hubbard model on a square lattice 

H = 'y ) tij a\,^aja + — 'y ^ NirjNia + Hc,ep, (1) 

i 


where a\^ and Oia are the Fermi creation and annihilation 
operators for electrons with spin a 12 (a = ±1, a = —a) 
on the lattice site i, and Ni = a a 

the number operator, tij is the electron hopping parame¬ 
ters (the nearest-neighbor hopping parameter t = 0.4 eV 
is used as the energy unit). The on-site Coulomb interac¬ 
tion (Cl) is U. The intersite Cl Vij and electron-phonon 
interaction (EPI) gij are defined by the Hamiltonian: 

Hc,ep = VijNiNj+'^ gijN^Uj, ( 2 ) 

where Uj are atomic displacements in particular phonon 
modes. In the strong correlation limit, U t, the 
projected electron operators referring to the single and 
double occupied subbands, the HOs, should be intro¬ 
duced [l^ : 

at = aL(l - + aim,, ^ Xf + Xf ■ 

In terms of the HOs the model H]) reads 

H = +^2X+ E E" 

i,a i iAJcr 

+ + (T(Af Af+ H.c.)} + He,ep, (3) 

where ei = —p, is the single-particle energy, £2 = U — 2 g 
is the two-particle energy, and /r is the chemical poten¬ 
tial. The HO A“^ = \ia){ifJ\ describes transition from 
the state |i,/3) to the state \i,a) on the lattice site i 
where (a, /3) refer to four possible states: an empty state 
(a, /3 = 0), a singly occupied state (a, (3 = a), and a dou¬ 
bly occupied state (q;,/3 = 2). The number operator and 
the spin operators in terms of the HOs are defined as 

m = ^Ar + 2Af, (4) 

(7 

= xr, St = {a/2) [xr - xn (S) 
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The HOs obey the completeness relation X™ + 
So- which rigorously preserves the con¬ 

straint that at any lattice site i only one quantum state 
a can be occupied. From the multiplication rule for the 
HOs = SjSjX^^ follows their commutation re¬ 

lations 




S^j{Sf,^xf±Ssc.xf) , 


( 6 ) 


with the upper sign for the Fermi-type operators (such 
as X^'^) and the lower sign for the Bose-type operators 
(such as the number ([4]) or spin dS]) operators). 

The unconventional commutation relations (jS]) for HOs 
result in the so-called kinematical interaction introduced 
by Dyson in a general theory of spin-wave interactions 
[ll|. To demonstrate the role of the kinematical inter¬ 
action in the model let us consider an equation of 
motion for the HO X['^ = al^aiattig ■ 

ij^Xf = [Xf,H] = {U-^l)Xf + [Xf,H,^,p] 

+ ^ tu {B^L>xf^ -aB^l,Xr') 

I,a' 

- J2tuX°^Xr° + aXf^), (7) 

I 

where the Bose-type operators are introduced 

= {Xf^ + xr)5.,, + xr5,,^ (8) 

= {Xi/2 + a S^) 5a'(7 + S'^ 6a'a, 

Bf^,, = iN,/2 + aS^)6a'a-S:Sa'^, (9) 


We see that the hopping amplitudes depend on number 
and spin operators caused by the kinematical interaction. 
In phenomenological spin-fermion models a dynamical 
coupling of electrons with spin fluctuations is specified 
by fitting parameters (see, e.g., [1-Q), while in Eq. ([7]) 
the interaction is determined by the hopping energy tij 
fixed by the electronic dispersion. 


GENERAL FORMULATION 

We consider superconducting pairing in the Hubbard 
model in the hole doping region. In this case the 
chemical potential /i is situated in the two-hole upper 
Hubbard subband and is determined by the equation for 
the average number of holes, n = 1 + 5 = (Ni) > 1. 

To study the electronic spectrum and superconduc¬ 
tivity in the model we introduce the two-time anticom¬ 
mutator Green function (GF) [l^ expressed in terms of 
the four-component Nambu operators, Xia and xj^ = 
{Xf'^ Af° Af2 xy) for two subbands: 

G,,a{t-t') = {{X,ait)\Xlit'))). ( 10 ) 


To calculate the GF m we use the equation of motion 
method by differentiating the GF with respect to time 
t and t'. Using the projection operator method 0 we 
derive the Dyson equation for the GF (uni) 0: 

G<,(k,a;) = [cc7o-Ea(k)-QZ<,(k,c.)]-'Q, (11) 

where tq is the 4x4 unit matrix and Q = {{Xia,Xl^}). 
The electron excitation spectrum in the generalized 
mean-field approximation (GMFA) is determined by the 
time-independent matrix of correlation functions: 

E-W = (12) 

j 

The self-energy operator is given by the multiparticle GF, 

Qz.(k,..) = {{zt^ I (13) 

The irreducible operators = [Xia,H] — J2i ^UaXia 
is determined by the equation {{Z-^\Xj^}} = 0. 

To calculate the self-energy matrix m we use the 
self-consistent Born approximation (SGBA) for the corre¬ 
sponding time-dependent multiparticle correlation func¬ 
tions. Assuming an independent propagation of Fermi- 
type excitations X[ ^, and Bose-type excitations Biaa' 
on different lattice sites we write the time-dependent mul¬ 
tiparticle correlation functions as a product of fermionic 
and bosonic correlation functions: 

B^j^^,\Biaa'{t)Xi 

= \B,aa' (t)) (14) 

The time-dependent single-particle correlation functions 
are calculated self-consistently using the corresponding 
GFs. This approximation results in a self-consistent sys¬ 
tem of equations for the self-energy (fOl) and the GF (fTTl) . 

The GFs for two subbands 1(2) in the normal state in 
the imaginary frequency representation can be written as 

{Gi( 2) (k, UJn)}~^ = iuJn “ £1(2) (k) - E(k, UJn), (15) 

where £ 1 ( 2 ) (k) are the quasiparticle energy (fT^ in 
GMFA. The self-energy for the two subbands can be ap¬ 
proximated by the same function: 

S(k,a;„) = -■^y]y]A(+)(q,k-q I w„-a;™) 

q m 

X [Gi(q, Wm)-f G 2 (q,a;m)] 

= iuJn[l- Z]^{uJn)]+ Xk{uJn). (16) 

In Fig. 1 we show the doping dependence of Z(q) = Zq{0) 
which weakly depends on 6 in the underdoped case for 
6 0.15 but sharply decreases in the overdoped region. 
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FIG. 1: (Color online) Doping dependence of the renor¬ 
malization parameter Z{q} along the symmetry directions 
r(0,0) ^ M{n,TT) X{n,0) r(0, 0) at T « 140 K 

for S = 0.05 (red solid line), S = 0.10 (blue dashed line), 
S = 0.15 (pink squares), S = 0.25 (black dash-dotted line), 
and 5 — 0.35 (black diamonds) [l^. 



FIG. 2: (Color online) Tc{5) in the WCA induced by all in¬ 
teractions (red solid line) and only by the spin-fluctuation 
contribution Xs/(blue dashed line) or only by the EPI 14p 
(black dash-dotted line) di- 


GAP EQUATION AND Tc 


Using the equation for the anomalous (pair) GF we de¬ 
rive equation for the superconducting gap function [14 1. 
In the linear approximation, for the gap in the two-hole 
subband v^k(w) = cr</52_cr(k, w) we obtain the equation 

[1 - 6(q)]^y(q,a;^) _ 

N ^ [uJmZqiuJmW + [£2(q) + Xq{uJrnW 




.)=§i: 


c{J(k-q) - U(k - q) + ^(q, k - q | - w^)}, (17) 


where 6(q) takes into account the hybridization effects 
of the two Hubbard subbands. J(q) and U(q) are the 
exchange interaction and the inertsite CL The frequency- 
dependent interaction is given by the function 

A(=^)(q,k|i/„) = -\t{q)\'^Xsf{K>^n) 
T[|U(k)|2xc/(k,:2„) -h XphiKi'n)], (18) 


where t(q) is the Fourier component of the hopping pa¬ 
rameter tij. The spectral density of spin (s/), charge 
(c/) fluctuations and phonons (ph) are given in terms 
of the dynamical susceptibility by the relations 0: 

Xs/(q,a;) = -((Sq|S_q)),^ , Xc/(q, w) = -{{6Nq\5N_q))^, 

and Xpft(q,w) = -((uq|u_q))<^. 

For comparison of various contributions to the pair¬ 
ing we consider the gap equation close to the Fermi 
energy, i^(k) = </j(k,w = 0). In this case instead 
of the dynamical susceptibility the static susceptibility 
x(q) = Rex(q, w = 0) appears in the gap equation; 




»’('>) tanhdSl 

2Ta 


[Z(q)]2 2E(q) 
x{j(k - q) - U(k - q) -t- |U(k - q)pXc/(k - q) 

+ l5(k - q)P Xph(M - q) ^(‘^o - |e(q)|) 

-|i(q)P Xs/(k - q)6»(ws - l^q)!)} , (19) 


where e(q) = £ 2 (q)/^(q) is the renormalized energy. 
Here ojq = O.lt and ojg = J = 0.41 are the cutoff en¬ 
ergies for phonon and spin-fluctuation excitations. Note, 


that in cuprates J(q) < U(q) and therefore the AF ex¬ 
change interaction J(q) cannot provide superconducting 
pairing proposed by Anderson [^. 

To solve Eq. (|TO)) we should introduce models for static 
susceptibilities. The spin susceptibility is determined by 
the function 


l+e[l + 7(k)]’ 


( 20 ) 


where 7 (k) = (l/2)(cosfca; -|- cosky) and ^ is the AF cor¬ 
relation length. The strength of the spin-fluctuation in¬ 
teraction is given by the susceptibility at the AF wave 
vector Q = (tt, tt), xq = Xs/(Q)j which is fixed by the 
normalization condition (S^) = (3/4)(l — d): 

^ 3(1 - _ 

2.., j7V^l-)-e^[l+7(q)] 



For the EPI we adopt a model with strong forward scat¬ 
tering proposed in Ref. [l5| . The static interaction in the 
model is determined by function: 

ffep(k) = |g(k)|2 xpft(k) = Pep (22) 


where the doping dependent parameter = l/(2i5) de¬ 
termines the radius of a “correlation hole”. 

To estimate various contributions in the gap equa¬ 
tion (IT^ we consider the d-wave model gap, ;p(k) = 
(A/2) {coskx — cosky). At hrst we consider solution of 
the gap equation (IT^l) for in the weak coupling ap¬ 
proximation (WCA), Z(q) = 1. As shown in Eig. 2, 
the largest contribution comes from the spin-fluctuation 
pairing induced by the kinematical interaction which is 
given by the averaged over the Fermi surface constant 
9s} = (kq)P Xs/(k - q))FS ~ 4t - 2 eV. Charge fluc¬ 
tuations and EPI contributions appear to be small since 
only the angular momentum Z = 2 of these interactions 
give contributions to the d-wave pairing 141. 
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FIG. 3: (Color online) Tc{S) for V = 0.0 (bold red line), 
V — 0.5 (blue dashed line), V — 1.0 (black dash-dotted line), 
and V = 2.0 (green dotted line) for t/ = 8 p^ . 


In the strong coupling approximation (SCA) in 
Eq. (ITTl) values of Tc are reduced by an order of mag¬ 
nitude as shown in Fig. 3 in comparison with the WCA 
in Fig. 2 due to large values of renormalization parameter 
Z(q)) in Fig. 1. Tc dependence on the nearest neighbor 
intersite Cl V in Fig. 3 reveals that the d-wave pairing 
survives as long as the Coulomb repulsion V does not 
exceed the kinematical interaction of the order of the ki¬ 
netic energy, V < it, 



FIG. 4: (Color online) Isotope exponent a{5) for g^p = 5t 
(red solid line) and g^p = 2.51 (blue dashed line). 

To study the isotope effect on Tc within the gap 
equation m we consider the mass-dependent phonon 
frequency ujq = y' Mq / {Mq + AM) = a;Q°^(l — /3). 

We neglect the polaronic effect for the d-wave electron- 
phonon coupling constant (|22l) (for discussion see 
Ref. [13 )• The result of numerical solution of the 
gap equation (HU for the oxygen isotope exponent 
(/? = 1/16) a = —d log Tc/d log M is shown in Fig. 4 
for two values of EPI gcp = 5t and gep = 2.51 in 
Eq. ((22ll . In accordance with an analytical estimation, 
a = {1/2) gcp/{gep + dsf), cr increases for larger values 


of gep ■ The doping dependence of the exponent agrees 
with experiments (see, e.g., 18|): it is quite small, 

a = 0.09 — 0.18, close to the optimal doping while 
drastically increases in the underdoped case at 5 < 0.1, 
a = 0.38 — 0.68 for gep = 2.5t — 5t, respectively. Similar 
results were obtained in Ref. 0) within the t-J model 
with EPI. 


To summarize, we have shown that in the limit of 
strong correlations a new coupling parameter, of the or¬ 
der of the kinetic energy of electrons, appears in the two- 
subband regime for the Hubbard model and brings about 
d-wave superconductivity with high-Tc- This kinemati¬ 
cal interaction is lost in the spin-fermion models and the 
slave-boson (-fermion) models treated in MFA. 


* Electronic address: plakida@theor.jinr.ru 

[1] N. M. Plakida, High-Temperature Cuprate Superconduc¬ 
tors (Springer, Berlin, 2010), 570 pp. 

[2] P. W. Anderson, Science 235 , 1196 (1987); P. W. An¬ 
derson, The theory of superconductivity in the high-Tc 
cuprates (Princeton University Press, Princeton, 1997). 

[3] M. Le Tacon,et ah. Nature Phys. 7 , 725 (2011). 

[4] M. Le Tacon, et al., Phys. Rev. B 88, 020501 (2013). 

[5] P. Monthoux and D. Pines, Phys. Rev. B 49, 4261 (1994). 

[6] T. Moriya and K. Ueda, Adv. in Physics 49, 555 (2000); 
Rep. Prog. Phys. 66, 1299 (2003). 

[7] A. V. Chubukov, D. Pines, and J. Schmalian, in: 
The Physics of Conventional and Unconventional Su- 
perconduetors, Eds. K. H. Bennemann and J. B. Ket- 
terson (Springer-Verlag, Berlin, 2004), Vol. I, p. 495; 
Ar. Abanov, A.V. Chubukov, and J. Schmalian, Ad¬ 
vances in Phys. 52 , 119 (2003). 

[8] Ar. Abanov, A.V. Chubukov, and M.R. Norman, Phys. 
Rev. B 78 , 220507(R) (2008). 

[9] J. Hubbard, Proc. Roy. Soc. (London) 276, 238 (1963). 

[10] J. Hubbard, Proc. Roy. Soc. (London) 285, 542 (1965). 

[11] F. Dyson, Phys. Rev. 102 , 1217 (1956). 

[12] D. N. Zubarev, Sov. Phys. Uspekhi 3, 320 (1960). 

[13] N. M. Plakida, in: Theoretical Methods for Strongly Cor¬ 
related Systems, Eds. A. Avella and E. Mancini (Springer, 
Berlin, 2011) Chap. 6, pp. 173-202. 

[14] N. M. Plakida and V. S. Oudovenko, Eur. Phys. J. B 86, 
115 (2013). 

[15] R. Zeyher and M. L. Kulic, Phys. Rev. B 53, 2850 (1996). 

[16] N. M. Plakida and V. S. Oudovenko, JETP 146, 631 
(2014). 

[17] N. M. Plakida, Phys. Scr. 83, 038303 (2011). 

[18] R. Khasanov, A. Shengelaya, E. Morenzoni, et ah, 
J. Phys.: Cond. Matter 16, S4439 (2004). 

[19] E.I. Shneyder, and S. G. Ovchinnikov, Pis’ma J. Exp. 
Theor. Fiz. 136 , 1177 (2009). 







